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We obtained an unprecedentedly large number of s-wave neutron widths through 7?,-matrix analy- 
sis of neutron cross-section measurements on enriched Pt samples. Careful analysis of these data re- 
jects the validity of the Porter- Thomas distribution with a statistical significance of at least 99.997%. 

PACS numbers: 24.30.Gd, 24.60.Dr, 24.60.Lz, 25.40.Lw 



O 



(N 



X 



> 

i> 

\o 
en 

o 
o 



% 



Neutron resonance parameters remain some of the 
most important information for testing random matrix 
theory (RMT) [l|, even more than fifty years after such 
data served as the original impetus for its creation. To- 
day, RMT pervades the physics of virtually all complex 
systems and, in the nuclear physics arena, is most often 
invoked in studies of quantum chaos [2| . Given the con- 
served symmetries involved, RMT for the Gaussian or- 
thogonal ensemble (GOE) of matrices is expected to cor- 
rectly describe fluctuation properties of nuclear levels at 
relatively high excitation such as near the neutron thresh- 
old. It implicitly assumes that reduced neutron widths 
^An ^^ s-wave resonances A follow a Porter-Thomas dis- 
tribution (PTD) [3[, which had been anticipated before 
RMT emerged. Currently, the overwhelming consensus 
is that P°jj data agree with the PTD. However, there are 
problems with both the data and analysis techniques used 
in reportedly the best test of the PTD to date 4] that call 
these results into question [sll . Measurement and analysis 
techniques have improved considerably since then, so it 
is worthwhile to perform new tests of the PTD. 

In this Letter, we show that rich T^^ data extracted 
from high resolution neutron total and capture cross sec- 
tions of i92,i94p^ measured at the Oak Ridge Electron 
Linear Accelerator (ORELA) facility display a significant 
departure from the PTD. To our knowledge, this result 
represents the most stringent test of the PTD to date, 
and the observed disagreement could have far-reaching 
consequences. 

Loosely stated, the PTD is based on the assumptions 
that s-wave neutron scattering is a single-channel pro- 
cess, the widths are statistical, and time-reversal invari- 
ance holds; hence, an observed departure from the PTD 
implies that one or more of these assumptions is violated, 
and so could be very interesting. 

To make reliable conclusions regarding the validity of 
the PTD, it is important that the data set be as pure, 
complete, and large as possible. Perennial problems with 
neutron resonance data have been (i) contamination of s- 
by p-wave resonances and/or resonances of neighboring 
isotopes, (ii) obtaining enough resonances with known 
spin J to perform statistically meaningful tests, and (iii) 



missed resonances due to finite experimental threshold. 

In the present case, problem (i) was minimized because 
Pt is near the peak of the s- and minimum of the p-wave 
neutron strength functions {Sq/Si ~ 10), and because 
we made high resolution cross-section measurements on 
natural Pt and four samples enriched in ^^^Pt, ^^''Pt, 
i95pt and I'^ept. 

Problem (ii) was addressed by combining data for two 
target nuclei ^^^Pt and -'^^^Pt, containing the largest num- 
ber of resonances (158 and 411, respectively), and by the 
fact that all s-wave resonances have spin J = 1/2 for 
these even-even nuclides. 

Finally, as described below, the novel approach of using 
an energy-dependent threshold in the analysis helps to 
solve all three problems. As a result, our i92,i94p-j- jg^^g^ 
are at least as pure, complete, and large as all previous 
r^j^ data which have been used for testing the PTD. 

Details of the measurements can be found in Ref. [6| . 
The ORELA was operated at a pulse rate of 525 Hz, 
a pulse width of 8 ns, and a power of 7-8 kW. Capture 
measurements were made at a source-to-sample distance 
of 40.12 m with a pair of CgDg detectors using the pulse- 
height- weighting technique, and were normalized via the 
saturated 4.9-eV resonance in the ^^^Au(n,7) reaction. 
Total neutron cross sections were measured on a sepa- 
rate flight path via transmission using a ^Li-loaded glass 
scintillator at a source-to-detector distance of 79.83 m. 

The 7?.-matrix code SAMMY [7] was used to fit both our 
transmission and capture data and extract resonance pa- 
rameters. Resonance energies E\ and neutron widths 
^jFaii, where gj is the statistical factor for resonances 
with spin J and Txm were used in the subsequent anal- 
ysis described below. For even-even targets, s-wave res- 
onances have gj = gi/2 — 1, and hence, ^jPaii = ^x-a- 

An asymmetrical shape in the transmission data could 
be used to assign ^„ = resonances [6|, see Fig. [T] 
However, there remained many weak resonances, most of 
which are p wave, for which we could not unambiguously 
determine the -^n value. As shown below, the potential 
problem posed by these resonances has been surmounted. 

The need to use a threshold on observables ^jPaii, 
to guard against possible systematic errors due to in- 



strumental effects and p-wave contamination, was real- 
ized very early on [3| in using such data to test theory. 
Use of an £^„-dependent s-wave threshold of the form 

3/2 

Tq = a qEyi , where oq is a constant factor, is a key im- 
provement in our method compared to previous analyses 
(in which only energy-independent thresholds were used) 
for at least three reasons. First, p-wave contamination 
is eliminated equally effectively at all energies. Second, 
as shown in Fig. [I] our instrumental threshold T follows 
very closely this same energy dependence; thus, possible 
diffusiveness of the instrumental threshold can be sur- 
mounted equally effectively at all energies. Third, statis- 
tical precision of the analysis is maximized by allowing 
the largest number of s-wave resonances to be included. 

The PTD is a special case (degrees-of- freedom 1^ = 1) 
of the family of x^ distributions, the probability density 
function (PDF) of which is denoted hereafter as /(a;|zy). 
Because theory predicts fluctuations, but not average val- 
ues, it is necessary to scale the data to their average, or 
expectation, value; x -^ ^\n/^{^Xn]^ where E[»] denotes 
the expectation value operator. 



ergies E\. The expression for this PDF reads 



We began our fluctuation analysis by employing the 
maximum-likelihood (ML) method, which has been used 
since the earliest tests of the PTD, to estimate v and 
E[r°„]. We used threshold Tq shown in Fig. [1] (with 
factor flo given in Table HI , which was chosen to reduce 
p-wave contamination to very low levels. 

The joint PDF for statistical variables F^^^ and E\ is 
defined in a 2D region I given by inequalities E\ < E^i^^ 
and r°jj > To{Ex), where i?max is an upper limit of en- 
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FIG. 1: Energy-reduced values of observables ^jFaii for indi- 
vidual resonances of ^^^Pt. Values represented by full-circle 
points belong to resonances to which we assigned in = 0. For 
the meaning of thresholds To, Toi and T see the text. 
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The factor C, ensuring a unit norm of /i*^, is v- and 
E[F°j^]-dependent. The ML function was calculated from 
all no pairs [E^'^^ , F^'^^] , obtained from the experiment, 
which fall into the region I. Specifically, 



L (^, E[F°„]) = n '^^ KT^ Ta^ I ^' E[FL 



(2) 



A contour plot of this ML function in the form 



z(^,E[F°J) = 2^ [lnL,„a, - \nL {i.,E[TU)] ^ (3) 

is depicted in Fig. [5] Here, Lmax is the maximum of the 
ML function. Results of the ML analysis are listed in 
Table H 

If In L displays near its maximum a shape close to 
a paraboloid, a contour for a fixed value z = k will en- 
circle approximately the ka confidence region of the ML 
estimates E[F°jJ and D (referred to hereafter as Poxp)- In 



this way, we found that the difference 1 



equals ap- 



proximately 2.7 cr for I92,i94pt and 1.4 cr for ^^^pt. The 
results for i92,i94p^ indicate that the PTD is excluded 
with high confidence. The result for ^^^Pt is in agreement 
with ^92,i94p^^ |-,ut with reduced statistical significance 
owing to the smaller number of resonances observed for 
this nuclide. For this reason, we did not include ^^^Pt in 
the subsequent fluctuation analysis described below. 

In the spirit of classical statistics, parameters v and 
E[FAn] are not random variables. Consequently, the func- 
tion z refers in ideal conditions to a distribution of esti- 
mates of these parameters, not to the parameters them- 
selves. Further, as can be deduced from Fig. [21 the 
shape of function InL (i^,E[r5^j^]) differs strongly from 
a paraboloid, which is also the case for ^^^Pt. So, at this 
point it is premature to draw a reliable conclusion from 
the values of Pcxp for I92,i94pt. 
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FIG. 2: Plot of z{i', E[r'^n]) constructed from the "^Pt data. 



TABLE I: Results of ML-based fluctuation analysis of s-wave reduced neutron widths of 



192,194,196 



Pt resonances. 



Sample 


-C/max 


flo 


no 


i^exp flOl Mo 


^ R S 




(keV) 


(eV-5) 




(eV-4) (meV) 


(meV) 


192p^ 


4.98 


7.00x10"* 


153 


0-57t"oit - 5.91±0.71 


(2.68±0.47)xlO"S" 26.1±4.5 0.9970 


194p^ 


15.93 


2.25x10"^ 


161 


0.47to;l8 9.4x10"'' 14.9±1.5 


(8.84±1.09)xl0"® 25.4±3.1 0.9975 


196p^ 


15.99 


3.19x10"^ 


68 


0-60to;26 9.4xl0"9 39.7±6.4 


(14.8±1.6)xl0"^ 



"Due to a high instrumental threshold, fj.^ was deduced assuming 
that 192, 194, 196 p^ share a common true value of /ig/ii-^. 



To check the veracity of the ML results, we under- 
took additional analyses. First, for a given target and 
a fixed value E[r^jj] we drew from the distribution gov- 
erned by the PDF, h°{Ex,Tl^ \ i/=l, Ep^J), a random 
sample, consisting of no pairs [E\,Tl^]. Then, with the 
aid of the ML analysis, we obtained an estimate P for 
this sample. From a large number of such samples we 
constructed the empirical cumulative distribution func- 
tion (CDF) of estimates P and determined a probability 
p^ — Pip > Pcxp I E[r^j^]) that a simulated value v is 
higher than the corresponding experimental value Poxp 
listed in Table [J Given the value of E[r°jj], the prob- 
ability Pu represents the statistical significance at which 
the validity of the PTD can be rejected. Values p^ ob- 
tained are plotted in Fig. [3] for -'^^^'^^^Pt. They are very 
high, but vary considerably with E[r"jj]. Therefore, we 
undertook further analyses to impose limits on Ep^^^]. 

To understand how this was achieved, consider a set of 
no independent variables {gi\ distributed normally with 
zero mean and unit variance. For such a set, the statistics 



i=l 



and Z^ 



Y.^9rf 



(4) 



will be governed by the same normal distribution, and by 
a x^ distribution with v = uq degrees of freedom, respec- 
tively. Therefore, for a set of values {g^^^} deduced in 
an appropriate manner from experiment, it is straight- 
forward to use Z and Z'^ to calculate probabilities for 
rejecting the null hypothesis that this set is consistent 
with the mentioned normal distribution. 

To employ these statistics, variable r° [^ needs to be 
transformed to variable g obeying the considered normal 
distribution. This was accomplished in two steps. First, 
the marginal CDF for the above-threshold widths 



rS E„ 



H%{TlnK]) 



/.°(srO'|i,E[r°])di?drO', (5) 



was used to transform Fj} to the variable r, which follows 
a uniform distribution. This was achieved by substitu- 
tion r = i/po (Fn I E[F°]). In the second step, with the aid 
of the inverse CDF of the normal distribution with zero 
mean and unit variance, G~^{r), we made a transforma- 
tion r ^ g, specifically hy g = G~^{r). Then, quantities 
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FIG. 3: Values of probabilities P(P > Pcxp), P(^>^cxp) and 
P (Z^ < Zcxp) as functions of EfF'Jj,] deduced from data on 
neutron resonances of ^^^'^^*Pt. Values referring to statistics 
V, Z^ and Z are plotted by symbols •, v a-nd '^i respectively. 



r, = if{^o(F°f^P|E[FOj) for z = l,2,...no, and hence 

sets [gl^"^} were calculated using values F^.^^ > Tq{E\.) 
from the experiments. Following this procedure, we eas- 
ily determined probabilities pz = P (^>^cxp | E[F°jj]) 
and pz2 = P {Z^ < Z^^^ \ E[TlJ) . Here, for a fixed value 
E[r°n], values of Zexp and Z^^^ are given by Eqs. (4) af- 
ter replacement gi -^ g°^^- Hence, probabilities pz and 
Pz2 represent separate statistical significances for reject- 
ing the PTD at various values of E[F5j^]. Both probabil- 
ities, calculated from i92,i94p^ data, are plotted in Fig. 
[3l As seen, the critical values of E[F°jj] for testing the 
validity of the PTD range from 5.3 to 6.9 meV and from 
15.0 to 18.3 meV for ^^^pt and ^^-^Pt, respectively. 

Before proceeding further, we performed an additional 
analysis to verify that p-wave contamination was negli- 
gibly small. To this end, we performed ML calculations 
based on a hybrid PDF for a mixture of s and p observ- 
ables gjTxn obeying two separate PTDs: 



h°^{Ex, g.jTxn\l^Q, 1^1,0)-- 
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Here, /ig and fj,i stand for expectatfon values E[r°jj] 
and E [(7jrAn(leV/i5A)^/^] referring to s- and p-wave 
resonances, respectively, while 6 represents the J- 
independent resonance-density ratio pj- /pj+ ■ The PDF 
/i°^ is defined in the region limited from below by thresh- 

3/2 

old Tqi = flQi-En with opi ^ oq which is at the 
same time safely higher than experimental threshold 

3/2 

T — aE-n ■ We assumed that p^ is J- independent. 
For mass numbers A w 190 this is justified, as it holds 
with about 1% precision that 51/2/53/2 ~ /(3/2)//(l/2), 
where /(J) is the resonance-density spin factor |9|. Fac- 
tor D ensures the unit norm of /i*^^. 

Following a path analogous to that described above, 
we constructed ML function L {(Iq, pi) and arrived at es- 
timates fi^ and /i^ given in Table HI In case of ^^^Pt, /i^ 
has been determined indirectly, see Table HI The method 
of its determination is supported by the fact that ra- 
tios ^/a^i for i^^'i^^Pt are within their rms uncertainties 
equal each other. In Table HI values of the dimension- 
less, energy-independent quantity R = {1 eV ' ) ao//i^, 
which represent ratios of To{En) to local expectation val- 
ues EpAn] for p-wave resonances, are listed for i92,i94p-j. 
targets. From these values and their rms uncertainties, 
we calculated probabilities of 0.069% and 0.0047% for 
-•^^^Pt and ^^*Pt, respectively, that among the observ- 
ables gjTxn > To{E\) there occurs one which belongs to 
a p-wave resonance, thus verifying that p-wave contami- 
nation is negligible. 

With the question of p-wave contamination settled 
then, the data in Fig. [3] indicate that for s-wave reso- 
nances in ^^^Pt and ^^'^Pt, the vahdity of the PTD is 
rejected with statistical significance levels S'192 = 0.9970 
and S'194 — 0.9975, respectively, in excellent agreement 
with our initial ML analysis. Because results for the two 
isotopes should be independent, the combined probabil- 
ity that the PTD is vahd is less than 1.2 x 10"^. Al- 
though we have shown above that it is very unlikely, if 
a p-wave intruder occurs among the s-wave ^^^Pt ob- 
servables, we calculate that if the resonance having the 
smallest g,jT\n value (relative to threshold Tq) is consid- 
ered to be p wave, this probability will increase in the 
worst case to 2.8 x 10"^. 

We conclude that our data reject the validity of the 
PTD with a statistical significance of at least 99.997%. 
This inescapable conclusion has been made thanks to rich 
experimental data obtained using state-of-the-art neu- 
tron spectroscopy, and the implementation of a novel ap- 
proach for testing the PTD. On the other hand, equally 
convincing evidence that the PTD holds for some or the 
majority of heavy and intermediate-weight nuclei is still 
missing. 



This result implies that at least one of the three as- 
sumptions behind the PTD is violated. For energies of 
the measurements reported herein, only elastic scatter- 
ing is possible, so the single-channel assumption is valid. 
Also, addition of another channel would result in i/ > 1, 
which would disagree even more strongly with the data 
than the PTD does. Violation of time- reversal invariance 
also implies v > 1, and therefore also is excluded. Hence, 
our results indicate the assumption that the widths are 
statistical is violated. However, there are no indications 
of non-statistical effects such as doorway states in the 
data. 

One possible explanation is suggested by the calcula- 
tions of Ref. [lO[ in which it was found that transition 
strength distributions deviated further from the PTD, in 
the direction of smaller j/, as model quantum-mechanical 
systems became more collective. Hence, our result that 
V « 0.5 for the ^^^^^^^Pt-|-n systems suggests the surpris- 
ing conclusion that i93^i9Spt display regular, rather than 
the expected chaotic, behavior at relatively high excita- 
tions near the neutron threshold. 

Alternatively, our results also could be interpreted as 
indicating that the PDF for reduced neutron width am- 
plitudes for the ^^^'-'^^'^Pt+n systems are not form invari- 
ant. In Ref. Ill| it was shown that this form- invariance 



assumption could replace the original [;3| somewhat qual- 
itative " statistical" assumption as part of a more general 
derivation of the PTD. Violation of this assumption could 
have far-reaching consequences. 
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